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Abstract 

The DD* scattering amplitude and the production of the final states tt + tt~ J /ip and 
7T + 7r~7r° J/ip near the D°D*° threshold are discussed following the recent suggestion 
that the observed peaks X(3872) and X(3875) in the decays B — > X K are due to 
a virtual state X in the D°D*° channel. The strong interaction is treated using the 
small interaction radius approximation. It is shown that the mass difference between 
the charged and neutral Z)W mesons results in a distinctive behavior of the relevant 
isotopic amplitudes. In particular, the shape of the peak in the tt + tt~J/^ channel 
should be significantly narrower than in the 7r + 7r~7r° J/'0 channel, which property can 
be used for an experimental test of the virtual state hypothesis. 



The narrow peak X(3872) [TJ [2j [3j H] currently commands a great interest due to the 
suggested possibility that it is dominantly a molecular state [SJ El [7] of charmed mesons 
(D°D*° + D*°D°), probably of the type discussed in the literature long time ago [El IS]- Such 
interpretation, as opposed to considering X(3872) as a regular charmonium resonance, is 
essentially based on two remarkable observations: the extreme proximity of the mass of X to 
the D°D*° threshold and the apparently strong violation of the isotopic symmetry indicated 
by the co-existence of the decays X(3872) -> tt+tt" J /iff and X(3872) -> 7i + 7i-7i°J/ip PEUD]. 
The measurement of the position of the X relative to the D°D*° threshold has been recently 
improved by the CLEO result [12] for the D° mass, which places the D°D*° threshold at 
3871.81 ± 0.36 MeV, and corresponds to M(X) - M(D°D*°) = -0.6 ± 0.6 MeV. Due to the 
exceptional closeness of the peak to the meson-pair threshold one can expect that the mass 
difference between the pairs of charged and neutral mesons, A = M(D + D*~) — M(D°D*°) « 
8.1 MeV, should naturally give rise to a significant isospin breaking in the properties of the 
state X. 

Furthermore, recent experimental study of the B meson decays B — > D D n K [T21 HI] 
and B — > D°D 0, -f K [14] revealed that the invariant mass recoiling against the Kaon displays 
a significant enhancement with a maximum at approximately 3875 MeV, which is only about 
3 MeV above the D°D*° threshold. The observed events can all be in fact attributed to the 
process B — > (D°D*° + DD*°) K since no distinction between the D*° mesons and their decay 
products was done. Moreover, the yield of the heavy meson pairs within the above-threshold 
peak is about ten times larger that that of the tt + tt^J/iIj and 7r + 7r _ 7r J/^ channels at the 
peak of V(3872). It has been most recently argued [T5] that a very plausible explanation 
of the observed enhancement of the D°D*° production combined with the smaller observed 
V(3872) peak in the 7r + 7r~ J ftp channel is that both these phenomena are due to a virtual 
state [TH [T7] in the D°D*° channel. In this picture the observed peak in the n + ir~J/ijj and 
7i + 7i~ 71° J/ ijj mass spectra is in fact a cusp with a sharp maximum at the D°D*° threshold. 

In this paper such possible virtual state near the DD* threshold is discussed within 
the approximation of small interaction radius, and a consideration is given to the isospin 
properties of the scattering and production amplitudes in this energy region. This approach is 
similar to the previously pursued [18] 'universal scattering length' approximation, and differs 
in including the effect of the nearby threshold for charged mesons D + D*~. An interesting 
energy-dependent behavior of the isotopic properties arises from the mere fact of the mass 
splitting A between the two isospin-related and coupled DD* channels. In particular it will 
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be argued that the expected pattern of the isospin breaking is consistent with the observed 
relative yield of 7i + 7i^J/ip and 7T + 7T~TT J/ip at the peak. Moreover, it will be shown that 
the production amplitude for the 1 = 1 state n + iT~J/ifj in the considered approximation 
necessarily has a zero between the D°D*° and D + D*~ thresholds, thus reducing the apparent 
width of the cusp and putting it in line with the experimental limit [1J T < 2.3 MeV on the 
width of the peak in this particular channel. 

It is assumed throughout this paper that the quantum numbers associated with the 
discussed peaks X(3872) and X(3875) are J PC = 1 ++ , corresponding to the S wave motion 
in the C even state D°D*° + DD*° and the coupled channel with charged mesons D + D*~ + 
D~D* + . These two channels are referred, for brevity, as respectively n and c. Considering 
the nonrelativistic dynamics of the mesons it is convenient to place the origin of the cm. 
energy E at the threshold in the n channel: E = M(DD*) — M(D°) — M(D*°). The energy 
range of interest for the present discussion is from few MeV below the n threshold and up 
to the c threshold, i.e. up to E ~ A. In this range the scale of the cm. momentum 
(real and virtual) in either channel is set by -y/2/xA ~ 127 MeV, where \l ~ 970 MeV is 
the reduced mass for the meson pair. One can apply in this region of soft momenta the 
standard picture of the strong- interaction scattering (see e.g. in the textbook [19J), where 
the strong interaction is localized at distances r < ro such that ro y/2fiA can be considered 
as a small parameter. Some well known points of this time-tested approach are repeated here 
in order to adapt the same treatment to the situation with two closely spaced thresholds 
in isotopically related channels. Outside the region of the strong interaction the motion 
the n and c channel is described by the free S'-wave radial wave function. Considering for 
definiteness an energy value between the two thresholds, < E < A, one can write the 
corresponding wave functions (up to an overall normalization constant) as 

Xn(r) = sin(/c n r + 5), Xc{r) = £ exp(-/t c r) , (1) 

where k n = y/2jlE and k c = ^2fx(A — E), 5 is the elastic^ scattering phase in the n channel 
and the constant £, generally energy- dependent, describes the relative normalization and 
phase of the wave function for the two channels. 

The wave functions ([1]) should be matched at r w ro to the solution of the 'inner' problem, 
i.e. that in the region of the strong interaction. In the limit of small r all the complexity 

1 In this consideration the small inelasticity due to the ir + ir~ J/tjj and ir + ir~ir° J/ifj channels is neglected 
and will be included later. Also the small width of the D* mesons is entirely neglected throughout this 
paper. 
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of the 'inner' problem reduces to only two parameters. Namely, in the region of the strong 
interaction the n and c channels are not independent and get mixed. Due to the isotopic 
symmetry of the strong interaction the independent are the channels with definite isospin, 
1 = and 1 = 1, corresponding to the functions Xo = Xn + Xc and Xi = Xn — Xc and 
the matching parameters are the logarithmic derivatives — k and —K\ of these functions at 
r = tq. Using the assumption of small ro the matching condition for the functions from 
Eq-(0Q) can be shifted to r = 0, so that one can write the resulting matching equations as 

k n cos6-£K c k n cos5 + £n c 

sin + £ sin o — £ 

These equations determine both the scattering phase S and the constant £ as 

co a = _^ (3) 

with 

2« «l - « c «l - «c«0 
^eflf = j , (4) 



K — Ki 

£ = - sm5 . (5) 



and 



The nonrelativistic scattering amplitude in the n channel is therefore given by [19] 
and the scattering length a is thus found from the E = limit of this expression as 
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The whole approach considered here is applicable if the scattering length is large in the scale 
of strong interaction. A large positive value of a implies an existence of a shallow bound state, 
while a large negative a corresponds to the situation with a virtual state [19] . According to 
the estimates of Ref. [15] the required by the data scattering length in the problem considered 
is -(3 - 4) fm, corresponding to a negative and quite small indeed parameter k^(E = 0) ~ 
(50 - 60) MeV. 

The physical picture, consistent with a small K e ff, and which could be argued on general 
grounds [8] is that an attraction in the 1 = channel is strong enough to provide a small 
value of Kq, while the interaction in the 1 = 1 channel is either a weak attraction or, more 
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likely, a repulsion. In both cases the absolute value of K\ is large, i.e. of a normal strong 
interaction scale, with the sign being respectively negative or positive. Another, purely 
phenomenological, argument in favor of large \k\\ is that no peculiar near-threshold behavior 
is observed in the production of the 1 = 1 charged states, e.g. D°D*~. At large \ki\ the 
expression (jlj) simplifies and takes the approximate form 

K c g « 2k - K c . (8) 

Using this approximation, one can readily see that in order for K e g(E = 0) to be negative 
and small the parameter k q has to be positive and quite small: 



k < vW2 « 63 MeV. (9) 

It is interesting to note that in the suggested picture the interaction in the 1 = state 
is strong enough by itself to produce a shallow bound state in the limit of exact isospin 
symmetry, i.e. at A — > 0. In reality the isospin breaking by the mass difference between the 
charged and neutral charmed mesons turns out to be sufficiently significant to deform the 
bound state into a virtual one, i.e. to shift the pole of the scattering amplitude from the first 
sheet to the second sheet of the Riemann surface for the amplitude as a complex function of 
the energy E. 

One can also notice that within the approximation in Eq.(jHJ) the scattering amplitude 
can be written in the form 

f = ~^ — rX7T> ( 10 ) 

which corresponds to equal coupling of the virtual state to the n and c channels. Such 
behavior, assumed in Ref. [15] on the grounds of isotopic symmetry, turns out to be applicable, 
as long as \k\\ is large, even if the isospin breaking by the mass difference is essential. 

Thus far the presence of any inelastic channels was neglected in the discussion of the 
DD* scattering amplitude. Such channels certainly exist and include the observed ones 
it + it~ J /ip (pj/ip), TT + n^n°J/tp (uJ/if)), 'jJ/ip an d probably other, which are yet to be found 
in experiment. The inelasticity however appears to be reasonably small, as one can infer 
from the observed [TBI E] dominance of the D°D*° production in the threshold region, and 
can be parametrized by a small imaginary shift 27 of the denominator of the scattering 
amplitude in Eq.(JB]) or Eq. ([TUl) : 

1 1 

/= « • (11) 

K, e s + ik n + z 7 2k — k c + i k n + z 7 
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The latter expression is similar to the one used in Ref.[15] and differs in that a term in the 
denominator, quadratic in k n , being neglected, as appropriate in the considered here small 
interaction radius approximation. The previous analysis [15] , based on the Breit- Wigner 
type description, includes such quadratic term and concludes that its contribution is very 
small in the discussed region of parameters. Furthermore, the relatively small value of the 
inelasticity is estimated there in terms of the scattering length as Ima/Rea ~ 0.1 or less. 

The width parameter 7 in Eq.( fITi) is the total sum over the inelastic channels, and in 
what follows the cotributions of the uJ/tp and pJ/ip channels, 7^ and j p , will be addressed 
in some detail. It is further assumed here [201 [15] that the 'seed' decay B — > XK is a short- 
distance process, so that the yield in each final channel coupled to X is proportional to that 
channel's contribution to the unitary cut of the amplitude /. This implies in particular that 

B[B^(D°D*°+DD*°)K\ : B(B - ujJ/^ K) : B(B - p J/j> K) = k n \f\ 2 : lw \f\ 2 : 7p |/p 

(12) 

where the specific expression for \f\ 2 depends on the value of the energy E relative to the 
n and c thresholds, as given by Eq. ( [111) an its analytical continuation across the thresholds. 
Besides the energy dependence of the overall factor |/| 2 , the heavy meson channel contains 
the phase space factor k n , while for the uJ/ip and pJ/ip yields an additional dependence on 
the energy arises from the factors 7^ and j p . 

A certain variation of the width parameter 7^ for the 7r + 7r~7r° J / ip channel in the discussed 
range of energy is of a well known kinematical origin. Indeed, the central value of the mass of 
the uj resonance puts the threshold for the channel ujJ/ip at 3878.5 MeV, which corresponds 
to E w 6.7 MeV in our conventions, i.e. squarely between the n and c thresholds. Any 
production of the 7r + 7r~7r° J/ip states at smaller invariant mass is a sub-threshold process, 
possible due to the width T w of the uj resonance. In other words, the energy dependence of 
the width factor 7^ can be estimated as [201 [IS] 

7 . = |AJ 2 gS ) , (13) 

where is the amplitude factor for the coupling to the ujJ/ip channel and is the effective 
momentum of uj at the invariant mass M calculated as 

(«)/,^ f M - m J/^ m w r w dm 2 

Qes ( M )= m m )\l—2 2 F2 ( 14 ) 

Jm (m z - m^y + m^L^ it 

with the cm. momentum |<f(m)| found in the standard way: 

l9lHI = - 2M ' ( > 
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The lower limit tuq in the integral in Eq. lTHl) can be chosen anywhere sufficiently below 
m uj — r^, since the Breit-Wigner curve in the integrand rapidly falls off away from the 
resonance. 

Numerically, the effective momentum can be estimated as varying from approxi- 
mately 20 MeV to 50 MeV between the n and c thresholds, i.e. when E changes from E = 
to E = A. As will be argued here, the amplitude A w should vary only slowly in the considered 
energy range, so that the shape of the threshold cusp in the ujJ /ip channel is determined by 
the behavior of the scattering amplitude factor \f\ 2 and by the estimated kinematical effect. 

In the pJ /ip channel the expected energy behavior of the yield is quite different. If one 
writes the corresponding width factor 7 p similarly to Eq. (1131) as 

7p=I4>| 2 <$> (16) 

the effective momentum qf^ can be estimated a s varying only slightly due to the large 
width of the p resonance: qf^ ~ (125 -r- 135) MeV as the energy changes between E = and 
E — A. On the contrary, as will be argued, the amplitude A p should experience a significant 
variation in this energy range and in fact cross zero between the n and c thresholds. 

In order to argue the claimed properties of the amplitudes and A p it can be first 
noticed that the inelastic processes DD* — > ioJ/ip and DD* — > pJ /ip involve a rearrangement 
of the heavy and light quarks and therefore cannot be due to peripheral interactions at long 
distances, but are determined by the dynamics at a typical range of the strong interaction. In 
other words the amplitudes of these processes are sensitive to the behavior of the heavy meson 
pair wave function at short strong- interaction distances. Thus for the 1 = state ujJ /ip the 
amplitude is related to the short- distance part of the function \oi while the amplitude for 
the 1 = 1 channel pJ /ip is determined by the function xi- Assuming also as previously, that 
the B decays are also determined by short distances, one can write for the decay amplitudes 
the expressions 

A(B -> ujJ/ipK) = J F u xa{r)dr and A(B -> pJ/ipK) = J F pX i(r)dr , (17) 

where F^ and F p are the weight functions in the respective channels, and their support, as 
argued, is limited to typical strong-interaction distance^. In the limit of vanishing inter- 
action range these functions should each be replaced by a 5-function: F U)P — > $^ iP (5(r), so 

2 Certainly, at present any details of these functions are not known. It can only be mentioned that it is 
likely that F p i=s due to the p — to universality. 
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that, using the expressions in Eq.flTJ and the results in the equations d2J), (jlj) and ([H]) for the 
solution to the matching conditions, one finds 

A(B -> ojJ/ipK) = 2 Kl - Kc sin5 « 2 $ w sin5 

fCo ~l~ ^1 — 2 K c 



A(S -> pJ/V = 2 $ sin 5 ~ 2 $ p sin 5 , (li 



where the latter expressions for each amplitude are written in the limit of large One 
can readily see that in this limit the amplitude for production of the isoscalar state ujJ '/if> is 
finite, while that for the isovector state pJ/ip is suppressed by Furthermore, it should 

be noted that, generally, in the limit of small interaction range it would be incorrect to 
retain such suppressed term without taking also into account effects of finite range. In order 
to fix this deficiency in the second line in Eq. lfTBj) one should go beyond the r = point 
approximation for the second integral in Eq. lfTTj) and use the first two terms of the Taylor 
expansion of the function Xi( r ) at r = rather than only the first term: 

Xi(r) = sin5-£+(£; n cos5 + f k c ) r + 0(r 2 ) = 2 K ° ~ K ° (1- Kl r) sin 5 + 0(r 2 ) . (19) 

Then the improved estimate of the second integral in Eq.l fTTj) can be written in terms of the 
effective radius R of the weight function F p , 

„ f F D rdr , . 

R = — P £ , 20 
JF p dr K } 

as 

A{B -> pJ/i/iK) = 2 K ° - Kc (1 - KiR)<& p sin<5 « 2 (k - « c ) (— - R \ <$> p sin 5 , 

Kq -\- K\ — 2 K c \K\ J 

(21) 

where, as discussed, the parameter K\R should be considered as being of order one. 

Clearly, the factor sin 5 in each amplitude is proportional to the scattering amplitude 
factor /, which is accounted for separately in the yield of each of the discussed channels 
(Eq. (fl2l ). Therefore this factor should be omitted in the amplitudes and A p , and one can 
conclude that the amplitude A^ is a smooth slowly varying function of the energy inasmuch 
as K\ is large. Furthermore, the ratio of the pJ/ip and LoJ/ij) production amplitudes is 
estimated as 

A >- K °- Kc <l- Kl R)^. (22) 



Aj Ki — K c $ 
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Since the situation where the X peak is a virtual state corresponds to a small positive kq 
satisfying the condition flU]) , the amplitude A p described by Eq. (1221) should necessarily change 
sign between the n threshold, where k c = ^/2/iA, and the c threshold, where k c = 0. It can 
be noticed that the existence of the zero of the amplitude A p results from both first terms 
of the expansion (fl9l) for the function Xi( r ) being proportional to kq — k c . 
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Figure 1: The expected shape (in arbitrary units) of the virtual state peak in the yield of 
7T + 7T~J/ifj (solid) and 7r + 7r~7r J/ip (dashed) channels. 



The expected difference in the shape of the cusp in the pJ/ip and uoJ/ip channels is 
illustrated in Fig.l. In these plots the parameters of the virtual state correspond to the 
scattering length a = — (4 + 0.5 i) fm, which is close to the possible fit values of the scattering 
length found in Ref . [15] . In the limit of large K\ this value of a translates into kq ~ 38 MeV 
and 7 ~ 6 MeV. One can see from Fig.l that due to the discussed zero of the amplitude the 
peak in the 7r + 7r~J/?/> is expected to be quite narrow in agreement with the experimental 
limit on the width of X(3872). The plots in Fig.l are normalized to the same total yield 
in each channel over the shown energy range in order to approximate the experimentally 
observed relative yield. Such normalization corresponds to setting 
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which value does not appear to be abnormal, even though at present we have no means of 
independently estimating this quantity. 

Clearly, the suggested significant difference of the shape of the virtual state peak in the 
two discussed channels provides with a way of discriminating between the options of X being 
a virtual state and a bound state. 

I acknowledge enlightening discussions of general properties of a low-energy scattering 
with Arkady Vainshtein and helpful communications with Yulia Kalashnikova and Alexander 
Kudryavtsev regarding their paper [Hj. This work is supported in part by the DOE grant 
DE-FG02-94ER40823. 
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